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ABSTRACT 


In this thesis we describe a pulsed NMR spectrometer 
suitable for work on solids. A novel feature is a simple 
tapped-coil arrangement that enables an impedance match 
to be obtained between the NMR tuned circuit and the 
equipment feeding it. Also we describe a sample holder 
that enables the orientation of single crystals within 


the NMR coil to be varied about two perpendicular axes. 


The spin-lattice relaxation time T, @£ oye ina 


doped saturated aqueous solution of NaNO, and in a single 


crystal of NaCl were measured in order to study the 


behaviour of the spectrometer. The Ty) of aa in the 


doped saturated aqueous solution of NaNO. was found to be 
20.4 ms, in rough agreement with a crude theoretical 


23 


estimatesantThe T. of Na “in: NaCl at. 299K was found to 


i 
be 14.620-15us.7 Thisa@issin onlystartreacreementswith 
other published values. 
Finally, the orientation dependence of the T) OL 


Se in a single, crystalo. NaNO, was measured. It was 
found to be independent of crystal orientation, as pre- 
dicted theoretically. Also, the T; value of 6.19 s at 299K 
is in excellent agreement with values found by other 
workers. However, such excellent agreement is fortuitous, 
since recent measurements with our equipment have shown 


that different values of T are obtained, depending on the 


control settings of the rf phase - sensitive detector. 
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CHAPTER 1 
INTRODUCTION 

Nuclei with spin number I 2 1 possess an electric quadru- 
pole moment in addition to a magnetic dipole moment. It 
has been shown by Pound (1950) that the interaction of the 
nuclear electric quadrupole moment with the time-averaged 
electric field gradient at the nuclear site affects the 
shape of the NMR absorption signal. In the case of a 
Single crystal, the interaction causes the NMR line to be 
split into several components. For a polycrystalline sample, 
it broadens the NMR line. Pound (1950) also showed that 
the interaction of the nuclear quadrupole moment with a 
time-dependent electric field gradient causes spin-lattice 
relaxation. “The main contribution to the Spin-lattice 
relaxation is caused by the fluctuating electric field 
gradient associated with phonons (Van Kranendonk, 1954; 
Wikner et al., 1960; Weber, 1963; Van Kranendonk and Walker, 
1968). However, the fluctuating electric field gradient 
associated with the rotation of molecular or ionic groups 
and with the diffusion of charged defects (Cohen and Reif, 
19577 Soaton, = 1965; Becker.” 1975) “through the lattice also 
contribute significantly to spin-lattice relaxation in 
certain temperature ranges. 

For nuclei such as 23N8 in NaCl which are situated 
at sites possessing cubic symmetry, the time-averaged 
electric field gradient is zero. The nuclear Zeeman levels 


are then equally spaced and only a single NMR line 
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is Observed. “Also, mutual spin flips or’spin exchanges of 
ereype i MY Mel «mt occur sulriclentlLy rapa) y to 
maintain a Boltzmann distribution among the nuclear Zeeman 
levels (Abragam and Proctor, 1958). In that case, the 
approach to equilibrium of Mo the component of the nuclear 
magnetization along the external magnetic field direction, 
Poeexponehntlat GMiener.- 1960." Zax, 1964)" i.e .*1s5° of the 


-t/T 


form M, =(M),C1-e : where T, is the so-called spin- 


1 
lattice relaxation time and (M >) (is the equilibrium value 

of M,- For nuclei” Situared at sites with’ lower thar cubic 
Symmetry, the spin-lattice relaxation is usually non- 
exponential. For example, for a spin system with I = 3/2, 
the, approach to equilibrium is in general given by the sum 
of 3 exponentials (Andrew and Tunstall, 1961). The three 
relaxation times are 1/2W, , 1/2W. and 1/2(W,+W.), where 

Wi and Wo are the probabilities of Am = +1 and +2 transi- 
tions respectively (Yosida and Moriya, 1956). 

Because of the non-exponential behaviour, experimental 
studies of the spin-lattice relaxation of nuclear spins in 
non-cubic single crystals are not easy to interpret. There 
fe ar so. the. adGgitlonal proolem that’ to 1S 7di PercuLt to 
selectively excite, by short rf pulses, individual NMR 
lines in multi-line spectra. Most measurements of the 
spin-lattice relaxation time in non-cubic single crystals 
have therefore been made with the crystal oriented relative 


to the external magnetic field so that the quadrupole 


splittings are essentially zero, i.e. so that all the 
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resonances overlap. When this is done, a Boltzmann 
distribution among the Zeeman energy levels is maintained 
by spin exchanges (Abragam and Proctor, 1958), just as in 
the case Ol a Gubicicrystal . es rhe spin-lattice relaxation 
behaviour is then exponential, and is characterized by a 
Single quadrupole relaxation time TT): 
The orientation of an external magnetic field rela- 
tive to a crystal has 2 degrees of freedom. The requirement 
of resonance overlap removes one of these degrees of 
freedom. The question can therefore be asked whether the 


quadrupole relaxation time T, is uniquely determined by 


a 
the overlap condition or whether it varies with the 
remaining degree of freedom. 

Bonera etial. (1970) and Avocadro .ct ale /C 1971). have 
measured the temperature dependence of the spin-lattice 
relaxation time of SNES in NaNO... Kasahara et ali. ClLon7) 
have measured the temperature dependence of the spin-lattice 
relaxation time of oon in NaNO.. None of these workers 
specified the orientation of their crystal with respect 
to the external magnetic field, except to say that. the 
erystal was set so that the resonances overlapped. | It 
would seem that these workers believed that TS is uniquely 
determined by the resonance overlap condition. 

Hughes and Spencer (1979) have shown that for certain 
types of relaxation processes, in particular the Raman- 


type two-phonon process, Ty) can be expressed as a linear 


combination of the components of a second-rank tensor, 
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subject of course to the constraint that resonance overlap 
occurs. Tt? can ‘then’ be shown ‘that T, is independent of 
crystal orientation for nuclei situated at sites possessing 
tetragonal, trigonal, hexagonal, and cubic point-group 
symmetries. This is not the case for nuclei situated at 
sites possessing triclinic, monoclinic and orthorhombic 
point-group symmetries. 

The 23Na iUCi es eT! NaNO. are situated at sites 
possessing trigonal symmetry. The spin-lattice relaxation 
time of these nuclei should therefore be uniquely deter- 
mined by the resonance overlap condition, provided that the 
relaxation is caused by a Raman-type two-phonon process. 
The aang nuclei is NaNO. are Situated at sites possessing 
orthorhombic symmetry. Thus, irrespective of the type of 
relaxation mechanism involved, the T) of these nuclei is 
not uniquely determined by the resonance overlap condition. 

The purpose of the work described in this thesis is 
to investigate the theoretical prediction of Hughes and 
23 


Svencer "(1979)" in *partrcurar=co “contirm *that T, for Na 


in NaNO,, is independent of crystal orientation. This spin 


3 
system was chosen in preference to the more interesting 


case of 23Na in NaNOo, because more complete resonance 


overlap can be achieved for oon a in NaNO, on account of the 
smaller second-order quadrupole shifts. In other words, 
the easier case was treated first. 


Spin-lattice relaxation times were measured by the 


so-called 180°-t-90° pulse method. In this method, a 
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"180° pulse" is applied to the nuclear spin system. This 
reverses the direction of the nuclear magnetization. A 
"90° pulse" is applied to the spin system at a time T 
later. This brings the nuclear magnetization into a 
transverse direction so that it precesses around the 
magnetic field and can then be detected as a rf voltage 
induced in the NMR coil. The amplitude of the NMR signal 
immediately following the 90° pulse is a measure of the 
value of MS immediately before the pulse. The spin-lattice 


relaxation time T, is obtained from the t dependence of 


1 
the amplitude of the NMR signal. 

The pulsed NMR spectrometer, which was built as 
part of this work, is described in the thesis. Also 
described is the method of data collection and analysis 
which avoids many types of systematic errors. Since this 
is the first time that pulsed NMR measurements have been 
made in this laboratory, some time was spent studying the 


benaviour, in. particular ene lLimrtbarions, of the spectro-— 


meter. These also are described in this thesis. 


a 


uine 


1 : may 
ty bawese ie see aed rane iin ox 


Hey aHtoein si) wp: gives eu ator aid 


ine cen Kgedinun oan Aen 


> 


mee 
4 9 ax befoseae cod mdi yong yn heh 


ch aly Ae . panei ” cy ma ok 
to sueseien .s.e5 wat ey? PO waft “jek ite ve a x Ma 
” ant - gelbo naly alo led fede keri am 
gal) 7 OFT mort Boeiiigse a io out? Ae iva at. 
to HIRE cald) De bard 


Bil 


’ 
5 
Leng & goidnstt: jab i Dedtan aay i 
Lee 
) s ‘he - 
sogte Sifaiateya. ty seta toe ab hom. tot 


4 atrswotzbeca ily besfig tela ont jg% + 


_ 
us goqn SAW ems QMoe , Yt bencid ait abat ¢ a 
se 7 Ss 
e Por 


} 
(4 %o Beiyas Pall etl yelonl ered er ' 


etdg @i baditese> otA chia. > She oy 


>= 


CHAPTER 2 


THEORY 


2.1 Description of NMR 


ook Energy: levelsdescuipulon 


If a nucleus possessing a magnetic moment ll and 
angular momentum nT is placed in a steady magnetic field H 
directed along the z axis, the so-called Zeeman Hamiltonian 


describing the interaction can be written as 


Yeo = Oe Ge) 


-yhl HO, C22) 


where tS is the z component of the nuclear spin i, and Yes 


the gyromagnetic ratio of the nucleus defined as 

y = u/al. (2.3) 
The eigenvalues of the Hamiltonian are 
EO = ~yfhimH |, (2.4) 
where m is the magnetic quantum number which may take any 


Ge tie, Zirh values i. let. 2. SL 


The energy levels of the nucleus are thus equally 
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spaced with the distance between adjacent ones being yfH . 
O 

Transitions between adjacent levels can be induced by an 

oscillating magnetic field perpendicular to Ho provided 


its frequency satisfies the equation 
ho. = yn C2). Oo) 
or 


U S= PYH 2. (2.6) 


The frequency Wo is called the NMR frequency. 

In thermal equilibrium, the lower energy levels are 
more highly populated than the upper ones, in accordance 
W2Ln the Boltzmann Gistribuvion. Since: the probability 
of an upward transition between a pair of levels is the 
same as that of a downward transition between the same pair 
of levels, there will be a net absorption of energy by the 
nuclear spin system from an oscillating magnetic field 
satisfying the resonance condition. This absorption of 
energy will tend to equalize the populations of the energy 
levels, a phenomenon called saturation. However, the 
interaction between the nuclear spins and other degrees of 
freedom of the sample material, the so-called spin lattice 
interaction, will tend to restore thermal equilibrium, If 
a non-equilibrium situation is produced, by a resonant 
oscillating magnetic field for example, the population 
differences will, in the absence of the oscillating 


magnetic field, approach their equilibrium values exponen- 
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tially with a time constant called the spin-lattice 
relaxation time TL: 

In practice, the Zeeman energy levels are not sharp, 
Since they are broadened by nuclear magnetic dipole- 
dipole interactions and other interactions. By measuring 
the rf energy absorbed by the nuclear spin system as a 
function of frequency, the so-called absorption line shape 
f(w) can be obtained. An expression for the mean square 
width or second moment of f(W) caused by dipole-dipole 
interactions has been given by Van Vleck (1948). 

In liquids, the dipole-dipole interactions are 
averaged out by molecular motions so that the absorption 
sienals are intrinsically very narrow. in such cases, 
the observed width and shape of the absorption signal is 


often governed by the inhomogeneity of the external magnetic 


fela iH, 
O 


2.1.2 Pulsed NMR 

An alternative description is the following. The 
elassical equation of motion for an object of magnetic 
moment i and angular momentum hi, situated in a magnetic 


=> 
field H, is 
6) 
du/dt = yu x H . Pr 


It can be shown (Slichter, 1963) that the expectation value 


=e . . 
of the magnetic moment <u> is given by a similar 
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equation, namely 
> > = 
d<u>/dt = y<u> x Ho ; G25) 


For a group of nuclei, the total magnetic moment or 


macroscopic magnetization is given by 
> > 
M = ) UL (229) 


where o is the magnetic moment of the kth nucleus. If the 
nuctel are identical; the equation of motion for the 


macroscopic magnetization is 
=5 > > 
GM/ dt. ="yMo x *H., G2 LO) 


It is convenient to view the motion of the magnetization 
from a frame of reference rotating about a. with an 
angular velocity Oi The equation of motion of the 
magnetization viewed from the rotating coordinate system 


is eCslichter, 1963) 


dM dt = M eal ne) zeal 
yas rot Y e) ). (2. ) 
If 3 = pei detay then dM ot / at is zero and Ca is time 


independent. This shows that the nuclear magnetization 
precesses indefinitely around H at a frequency 8 which 


is numerically equal to Wo (see Eq. (2.6)). The frequency 
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2 is called the Larmor frequency and is thus equal to the 
NMR frequency. The senseof the precession is governed 
by tne sign y. 


Suppose a rotating magnetic field 
H, = H,(icoswt + jsinwt) (202) 


directed at right angles to H is applied to the spin 
system. Viewed from a coordinate system that rotates 
about the z axis at the same frequency w as Hee and with 


Hy directed along the x axis, the equation of motion of 


the magnetization becomes 


dM. /dt = 


aes x[iyH, + k(wt+yH,) J. (25.13) 


M 
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In this rotating frame, the magnetization precesses about 


a static effective magnetic field 
> “A 
H an Os es k(H + —)); Cae 14.) 


At resonance, when w = ~YH,> i.e. when w is numerically 


equal to the NMR frequency Wo» 
= iH aa) 


and the effective field is directed along the x axis of 


the rotating coordinate system. If the nuclear 
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magnetization is initially parallel to Hi as it would be 
in thermal equilibrium for example, then the magnetization 
will precess in the yz plane’ about the x axis of the 


rotating frame with a frequency yH if the rotating 


L 


magnetic field H, is applied for a short time i say, 


1 
the angle through which M precesses in the rotating frame 
is given in radians by 


6 = yH C226) 


oe 
Piecwslpulacion as Shown Hnetig. 2.1 Ga). OL partieular 
importance are so-called 90° and 180° pulses whose lengths 
(or amplitudes) are chosen so as to turn the nuclear 
Magnetization through 90° and 180° respectively in the 
rotating coordinate system. The effect of such pulses 
arTeeshnown in. bacs. 2217 (b)7 and Gee 

Suppose a 180° pulse is applied to a spin system 
which was initially in thermal equilibrium. The pulse 
reverses the magnetization so that the z-component Me is 
equal to -M,> where My is the value of the magnetization 
in thermal equilibrium. Because of the spin-lattice 
interaction, the nuclear magnetization will relax toward 
Mo: For a spin system with I = 5 the relaxation is 
always exponential and the time constant is called the 
longitudinal relaxation time or spin-lattice relaxation 
time T,. The time dependent magnetization following a 


1 


180° pulse at time t = 0 is therefore 
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BigureeZont 
Precession of the nuclear magnetization (for y > 0) 
about the x axis of the rotating coordinate system when 
a rf magnetic field Hy is applved! alone the x axis ser 
a short time vas In case (a) yH.t =68> Cases 0b) sand 
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(c) show 90° and 180° pulses where ESS, = 71/2 and T 


respectively. 
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ae M[1-2exp(-t/T,)]. Cov 


Suppose a 90° pulse is applied to a spin system 
which was initially in thermal equilibrium. Immediately 
after the pulse, the longitudinal or z-component of the 
magnetization is zero and we have a transverse component 
of magnetization Mo which precesses around oe with an 
angular frequency Wo in the laboratory coordinate system. 
This precessing magnetization, with its associated time- 
dependent magnetic flux, induces a rf voltage in a coil 
wound aroundthe sample. This voltage is the so-called NMR 
Signal. Because of the magnetic dipole-dipole interaction 
between nuclear spins, the transverse magnetization decays 
toward its thermal equilibrium value of zero. The 90° 
pulse and the decaying NMR signal are shown schematically 
in Fig. 2.2 (a). The NMR signal, or more precisely its 
envelope indicated by the dashed line in the figure, is 
called the free induction. decay F(t) say. It can be 
shown (Abragam, 1961) that F(t) is the Fourter transform 
of the absorption line shape function f(w) mentioned in 
the previous section. In some cases, especially liquids 
and gases, this decay is exponential and the time constant 
is called the spin-spin relaxation time To. inisolads: 
the spin-spin relaxation is usually non-exponential. 
However, an effective To, is often used to indicate the 
time scale of the decay of the transverse magnetization. 


In liquids and gases, the decay is often dominated by the 
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Figure 2.2 
Schematic representation of (a) the FID following a 
90° pulse, (b) the FID following a 180°-tT-90° pulse 


sequence. 
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the inhomogeneity of the external magnetic field. 

The spin-lattice relaxation time T; can be measured 
by performing an inversion recovery experiment in which a 
180° pulse followed a time t later bya90° pulse is applied 
to the spin system. The 180° pulse reverses the z 
component of the magnetization. This then recovers toward 
Ptessequiiibrium vyaluesin, accordance: with’ Eq. (2.17). The 
z-component of magnetization, which is not directly 
observable by NMR, is made observable by applying a 90° 
pulse which rotates MS into the transverse direction. The 
pulse and the free induction decay are shown schematically 
in Fig. 2.2 (b). A measurement of the amplitude of the 
FID immediately following the 90° pulse is a measure of 
M. acla time t aiter the 180° pulse. By rearranging 


Z 


loge) oa.4), Lt can. be shown that 


gn(M-M,) = 2n 2M,-t/T, - C2413) 


Thus by plotting &n(M-M,) aoa tunCt ron OL tone r;cal 


determine T, from the slope of a straight line graph. 


1 


It can then be mooted that it 1S not’ essential for 
the magnetization to be perfectly reversed. If the 180° 
Hiulee 16) not, perlect, M, will still recover exponentially 


according to an equation of the form 


in(Mo-M_) =A - t/T,- 62-139) 


wall 
os rm 


i 


bertenws od ano an 


ee an | 


a doldw nt Jopme* 


7 
7 


Set coterie 
ie 
tot tat + omit a 


i Be - ) 


‘betiqgs et salaq * 


=» ona cooTaeat eating "OBI oat 


a 


pawns Ry afet feds erat wobdanktee om oft 
oan -_ 
ett CFE St) .@ del eoanb tooo: ut estat) 


- - Tee 


of, - 


, 


igondib shim #2 dotde forte sitanaaet we 
‘or »« apixiaan vd sidavisedo oban at ,4ae vd ‘w 
r eel ih eetovenett od? oft. en ‘dole 
* ip! 
i Peittod iwone wts g50pD nts? bart ett ons 
a¢€ weet cE 4 a7 1f : it jo tues Sem A . ba) S28 « 
: a 
enn « ef Safuq “0? ‘sd? yotiot lot. eleserae 
a 
‘woieot VA seedling “CORN ott eta? any Late 


secs nvede ed cas ws cnt Yo 


(ar.o) ep Tire ae = (> Oe 
io SeNoOR & ee (LM Me Bee 
,0kL Siyterse 6% oqgote ety sag y 

ie ‘ 
10. \eivaseap Om @ tt Jady beton od cat 
‘oar oft %] wvepdeves vildostivteq sf of oer 


el toecodRre teVaoot tizjdge tity _ 


r 


Gxbi a#i To solde 


Also, the 90° pulse need not be perfect, so long as 
ci T.5 which implies that T) ae To: In many solids this 


is the case, so that T) measurements should be free from 


Systematic errors caused by pulse imperfections. 


2.2 Nuclear quadrupole interactions 


FOr aaniclets wich nasea "spin number Ls 1) at is 
necessary to take into account the electrostatic inter- 
action of the nucleus with its surroundings. The electro- 


Static Hamiltonian is 
WK = | PCr av CaN (2520) 


where o(r) is the nuclear charge density, V(r) is the 
electrostatic potential of the charges external to the 
nucleus and the integral is over the nuclear volume. By 
expressing the potential V(r) in a Taylor expansion 


about the nuclear center of mass, the Hamiltonian becomes 
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> 1 
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(2.21) 


where | o(r)dv is Ze the nuclear charge, 


| p (Fx dV ie the nuclear electric dipole moment, 
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moment tensor. 

Since the first term in the Hamiltonian is orientation 

independent, it is omitted. Furthermore, because of 

parity considerations, the electric dipole moment (and 

the electric octapole moment etc.) relative to the 


nuclear center of mass vanishes. Thus, the expression 


for the Hamiltonian can be written as 
al, 
ww = 5) Q!,V., + hexadecapole term + 
ae JKOIK 


where V., = 9°V /ax 5 8x The hexadecapole interaction 


jk k 


is reeligible in all practical cases (Abragam,” 1961). 


Therefore we need only consider the quadrupole 


Hamiltonian 
1 
Gt : ; 2 
Ay 2 enue ( 
The quadrupole moment tensor is is symmetric and thus 
has six independent components. It is convenient to 


define a simple tensor ork which is both symmetric and 


traceless. This is done by writing 


= ‘1, = oy 
Q5% oni yr L Sees ( 
where 6., is the Kronecker delta function. By 


jk 
substituting Oak from Eoavet 2.24) an Eq. (2.23) and 


(2. 


av is Se the nuclear electric quadrupole 


22) 


,23) 


24) 
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ignoring the term involving traces since it is independent 


of the nuclear orientation, it is. found that 


Ha = e). Sa ee (2.25) 
By using the theorem that the corresponding matrix 
elements of all traceless, second-rank, symmetric tensors 
are proportional (Ramsey, 1953), it can be shown (Cohen 
and Reif, 1957) that the matrix elements of the quadrupole 
Hamiltonian are given by 


na es 2 
<m'|}{|m> = Leeueh Gat sa)] idem | 2 (151 ai 15)-6 51 |m>V. 


j kee jk 
C2526) 
where Q is a scalar quantity called the nuclear 
quadrupole moment and is defined as 
eQ= | o(%) (327-r7) av (2.27) 


where 0(Y) is the nuclear charge density and z is along 
the axis of the spin of the nucleus. By introducing the 


raising and lowering spin operators, 


Lepurl ast if (2.28) 


the Hamiltonian becomes 
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In a representation that diagonalizes the Zeeman energy, 


the matrix elements of Hy are 


<m|H|m> = Al 3m*-1(1+1) ]v, (2.31) 
<m+1 [FQ |m> = A( 2m#1) (zm) (I#m+1) ]? Vey C223) 
<m+ 2 [J | m> = AL( 1pm) (Igm-1)(I#m+1) (14m+2) ]? V9 C2933) 
<m' |, | m> =30° watox |m'-m|>2 (234) 


where A = eQ/4I(2I-1). 

In the so-called high field case, the quadrupole 
Hamiltonian is treated as a perturbation on the Zeeman 
energy. The Hamiltonian of the spin system is then 


given by 


gic 


oF 


utbeue ots ,2ab>o hfett Agia berkny-cm ” 
so oft #0 MORTAINMSG & one incnadlis 
ood! af sedeyve aiee dd) — 26 


a 


Tor eh pa] : 


2inmtele abated 
i 

asl 7) 

eh cret teas Se cant 
: > 


a i 


7 
‘Nit teats. = <0) glee 
5. aye 

oe 
a 


a tying rie } ~ ent gtleiae 


2 rt 


‘cs 


“ 


itr] 207 o od ie aie 

- 7 

(eats ONG ie . 
a 7 


ad Per 


a 7 - 


”v ag as 


Se gm ein | 2 5t, (2.35) 


where the matrix elements OTS are given in Eqs. (2.31), 
(Ze oele* G2a30) £4iee JO4) : 
The energy to the second order in the quadrupole 


interaction is (Bersohn, 1952) 
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The NMR frequencies are governed by the time-averaged 
values of En and these values are in turn governed by 

the time-average values of the components of the electric 
field gradient tensor. To take account of crystal 
symmetry, we choose a different coordinate system from 
that used previously. Let xyz be a principal coordinate 
system fixed with respect to the crystal so that 
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(2.37) 


It is also convenient to define two parameters q and n by 
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The electric field gradient tensor components obey 
Laplace's equation v-v =O.) -bnus, in Satisties the 
condition 

OA Seats le C2540) 


If the nuclei are situated on 3-fold, 4-fold or 
6-fold axes of symmetry, Vee and Ne are equal and n is 
zero. Therefore n is called the "asymmetry parameter". 

We denote the polar and azimuthal angles of the 
external magnetic field a Wiki respect. to, une, crystal 
coordinate system (xyz) by ®@ and $¢ respectively, as shown 
Peele. aio bY CLanStormane wha. 7 U2) 56 san. Omytne 
laboratory coordinate system to the crystal coordinate 


system, it can be shown that 
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Pyeure 2.0 
Orientation of the external magnetic field relative to 


the principal axes of the crystal. 
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where v, is the unperturbed frequency YH / 21. 
The frequencies of the quadrupole split NMR lines 
are given by 


eet en be (2.42) 
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For a spin system with I = 3/2, it can be shown by using 

Eqs. (2.41) and (2.42) that the frequency of the so-called 
centre® ine yreorresponding* torthe transition m*="1724<> — 1/2, 


is given by 
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Kiso the frequencies v~ of the so-called satellites, 


corresponding to the transitions m i= 3/2 «> 1/2 and 


mm =-= 1/2 +> = 3/2, are given by (for + positive y) 
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We see from Eqs. (2,43) and (2.44) that, to first order, 
the centre line frequency is not perturbed by the 
quadrupole interaction, and the satellites are equally 


spaced relative to the centre line. In this approximation, 
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it follows that the centre line and both satellites 


overlap when 


Sena Gal eacineoces26 = 0 C2545) 


However, it can also be shown from Eqs. (2.43) and (2.44) 
that, when account is taken of the second order terms, 
perfect resonance overlap will never occur. Thus, 
whether an approximate overlap of the centre line and 
satellites can be achieved, depends upon the relative 
magnitude of the second order quadrupole shift and the 


width of the individual NMR lines. 
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2.3 Orientation dependence of the nuclear quadrupole 
spin-lattice Tellaxation in Crystals. 


ExXamiInationeon Bis ee wis (2. oo). and Cero) 10 
the previous section shows that the nuclear quadrupole 
Hamiltonian couples Zeeman states whose m values differ 
ny ir andy +25 Using rerm. os Golden Rule,’ tie probabllicy 


of transitions where m changes by +1 and +2 is given by 
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In Eqs. (2.46) and (2.47), the electric field gradient 
Operators such as Mae are viewed from the laboratory 
coordinate system x'y'z' in which the external magnetic 
ftetd lies along the z' “direction. Also 6 is 1/kT, , TL is 
the lattice temperature, and n and n' are the initial and 
final lattice states respectively. It is convenient to 
define (cf. Yosida and Moriya, 1956) two quantities Wi and 


Wo, which characterize the nuclear quadrupole spin-lattice 


relaxation. They are defined by the expressions 


W, (2m+1)°(I+m+1) (I-m) 


W ee ae ee (2.48) 
te 21(2I1-1)" 
and 
Wo(I-m)(I-m-1)(I+m+1)(I+m+2) 
W SOT Greeny sae ic eee yep ge ae (2.49) 
tne STC 211) 
For the case I = 3/2) with which we are concerned in this 
thesis, ar is zero except when m equals -3/2 or 
Lia, oamilarly’, Ee Resid: is zero except when m equals -3/2 
Ore-1/2. The electric. duadrupole transitions for [= 3/2 


are shown in Fig, 2.45 “lt follows trom Eqs. (2.48) and 
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Figure 2.4 


Nuclear electric quadrupole transitions Oe oe 
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depend upon the direction of z', and hence the direction 
of the external magnetic field, relative to the crystal. 
It can be shown (Hughes, 1973) by means of the principle of 
time reversal invariance that the general orientation 


dependence of Wy and Wo is Of the torm 
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where 6 and o are the polar and azimuthal angles of the 


external magnetic field measured relative to a crystal 
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coordinate system. The coefficients eg are in general 
linear combinations of the components of a real symmetric 
fourth-rank tensor (Pietila, 1968; Snyder and Hughes, 


19 Tigetughes, 11973) 
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Qn a, ak n | a8 | | Gis) | ( n n mtu =m 
n 
) exp(-BE_) 
" (2.53) 
WHeLreo Op, ewe br = xX, Vi, 2. Mere xyZ- are the 


axes of a coordinate system fixed with respect to the 
crystal and chosen in such a way as to make full use of 
crystal symmetry. The tensor components Ma Bat grou em™ 
are properties of the crystal. In general there are 
fifteen independent M-tensor components. However, the 
presence of crystal symmetry reduces this number. For 
example, for 3 and 2m symmetry there are four and six 
independent M tensor components respectively. 

The approach to equilibrium of a spin system (with 
I > 1) is usually non-exponential. In general, the 
relaxation behaviour of a system of nuclei with spin 
number I is governed by 2I relaxation times (Andrew and 
Tunstaid. 19Giy.- Forcexampre, for I = 3/2, the approach to 


equilibrium is of the form 
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where A, B and C are determined by the initial conditions. 
However, if a Boltzmann distribution can be established 

among the Zeeman levels of the spin system, the relaxation 
behaviour is exponential and the nuclear quadrupole spin- 


lattice relaxation time T, is given by (Andrew and 


1 
Tunstall - 1961) 


ii = [(21+3) /51( 21-1) ](W,+4W,). (2954) 
Spin exchanges (Andrew and Tunstall, 1961) maintain 
a Boltzman distribution among the energy levels of a 
spin system if the various quadrupole-split resonances 
can be made to overlap. —As-—-pointed-—out-in section (22.2), 
this can be achieved for a system of identical nuclei in 
a non-cubic crystal by suitably orienting the crystal 
relative to the external magnetic field, provided second 
order effects are sufficiently small. 
For the case of relaxation mechanisms such as the 
two-phonon Raman process (Van Kranendonk, 1954; Van 
Kranendonk and Walker, 1967, 1968), the energy difference 


E ET is negligible in comparison with typical lattice 


m-u 
energies EF, im the: idelestartunction tine Fqsie(i2.60)6 <2.51) 


and (2.53). The M-tensor components associated with 


Wi are therefore essentially equal to those associated 


with Wo: Byiwsubstaitutineye ne Lois ¢250),, ta doertbhen: found 


(Hughes and Spencer, 1979) that 
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-l_e 29? (2I1+3) 
ah at vo: Ee [k, +k, 
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+sin“0(k,cos26+k,sin29)], C2555) 


where the coefficients k, are linear combinations of the 
M-tensor components. The angles 6 and @$ are not 
independent, since they must satisfy Eq. (2.45) for 
resonance overlap. By eliminating @ using Eq. (2.45) it 


can be shown that 


1/2 ; 
alee oO (2143) a +[2(1-ncos29)]°""(k,cos¢+k,sin$) 
; 801°( 21-1) ; 3 - ncos2¢ 
2(k-cos2ot+k,sin2¢) 
7 (2.56) 


3-ncos2¢ 


Because of crystal symmetry, some of the coefficients 


wants.” ‘nus, or one in single crystals of NaNO. where 
the point group symmetry is = it can be shown by using 


Ba. (2.06) that 


T.* = (e7Q7/60) (8 M ne: 


1 Tic M3131 + M3333) - (2.57) 


We note that T) is independent of the orientation of 


the crystal relative to the external magnetic field. This 


is not the case in general. For example for a in NaNO, 


where the point group symmetry is m2, T. is given by 
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CHAPTER 3 


EXPERIMENTAL DETAILS AND RESULTS 


own ADparabus 


bvbei) ‘Briet descriptionsorl the pulsed NMR spectrometer 


In this section we briefly describe the various 
parts of the apparatus and their functional relationships. 
Important parts of the apparatus are described in more 
Celatiw gn sections Sy tao to, Sil a7’ 

AuDLOeCK Gane raMm Ot tnereapparatus: as showns inekig. 3.1. 
A Takeda-Riken 3130 frequency synthesizer supplies a cw 
Signal to a Matec 515 rf gated amplifier. The gating 
Signal is provided by a pulse generator. The second 
connection between the pulse generator and the rf gated 
amplifier carries the reset pulses to be discussed later. 
The clock pulses for the pulse generator are derived from 
a cw Signal fed from the frequency synthesizer. This 
ensures that the gating pulses produced by the pulse 
generator are coherent with the rf signal being gated. 
The frequency synthesizer also provides a cw reference 
signal to the phase-sensitive detector. This signal is 
passed through a phase shifter in order to allow adjust- 
ment of the phase of the reference signal. The phase 
shifter uses a variable 10-turn delay line to provide the 
phase variation. 


The output from the rf gated amplifier feeds a 
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Block diagram of the system. 
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Single NMR coil in which the sample is located. A lumped 
circuit duplexer is used to decouple the NMR coil and the 
rf phase-sensitive detector when the rf pulses are on. 
When rf pulses are off, the NMR coil is coupled to the 

rf phase-sensitive detector so that the (weak) NMR signals 
can be amplified and detected. The phase-detected NMR 
Signal is fed to a Nicolet 1170 signal averager.* The 
averaged NMR signal can be viewed on a CRT, digitally 

read out via a Texas Instruments type 733 ASR electronic 
data terminal or read out on a Hewlett-Packard 7035B X-Y 


recorder. 


3.1.2 Sample probe 

The NMR sample coil has a length of 0.8 inches and 
a -dLameter ot OF 56" inches. It consists of"21 turns of 
Silver wire 0.0125 inches in diameter and has a measured 
inductance of approximately 7uH. Silver was chosen in 
preference to copper because the NMR frequency of copper 
is close to that of sodium, and it was important to ensure 
that no background NMR signals were present. 

The NMR coil, represented by the inductance Ly in 
Fic. 3.2, 1S connected to a variable 5-50 pF capacitor 
Ci: In order to avoid dielectric breakdown caused by the 
high voltages present during the rf pulses, a vacuum 
eapacitor from ITT Jennings Ltd. was used for C,- 
*I am grateful to the Department of Physics, University 


Of Gaivary. in particurar Dr, A.W. Harrison, for the: loan 
of the Nicolet 1170 signal averager. 
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The capacitor was physically small enough to be located 
adjacent to the NMR coil within the sample probe. This 
kept stray capacitance to a minimum and allowed the NMR 
coil to have as many as twenty-one turns. A ferromagnetic 
bearing on the capacitor was replaced by a nylon one so 
as not to disturb the homogeneity of the external magnetic 
field. 

The output impedance of the rf gated amplifier is 
902 whereas the impedance of the LjC, combination at 
nesonance is ‘several stens of thousands of ohms. An 
impedance match was obtained by connecting the rf gated 
amplifier across a single turn of the NMR coil as shown in 
Fig. 3.2. Connection across the single turn was achieved 
by soldering two lengths of silver wire 0.0125 inches in 
diameter to the coil. The tapped-coil behaves like a 
transformer, and the impedance ratio is thus approximately 
the square of the turns ratio, i.e. De or 441. The input 
impedance of the NMR coil (measured across the single turn) 
was found to be about 1002 at an operating frequency of 
14.5 MHz when C. was tuned to resonance. A resistor was 


lf 
therefore connected across Ly and C, and the value was 
chosen so that the input impedance was 502 and resistive 
at the operating frequency. The impedance measurements 
were made using a Hewlett Packard 4815A rf vector impedance 
meter.” Since the output impedance of the rf gated 
*I am grateful to the Department of Electrical Engineering, 


University of Alberta, for the use of the Hewlett Packard 
4815A rf vector impedance meter. 
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amplifier is also 502, the sample probe and gated 
amplifier are connected by ordinary 502 (RG-58) coaxial 
cable. This impedance match ensures maximum power 
transfer from the amplifier to the probe. 

To our knowledge, this is the first time such a 
tapped coil arrangement has been used. Most other NMR 
spectrometers use a matching network consisting of 
several components. The tapped coil arrangement is simple 
and is also capable of generating large rf magnetic fields 
with a relative modest power amplifier. 

it 1s important that che ‘@ factor of the NMR coil 
should not be too high. Otherwise, fast rise and fall 
times for the rf pulses cannot be obtained. The Q factor 
of the sample tuned circuit was found to be 34, and this 
is low enough for work on a el nuclei. 

The probe body was solidly constructed to reduce 
microphonic. effects. Brass was used in the construction, 
in preference to aluminum for example, to reduce spurious 
ringing signals observed by other workers (Buess and 
Petersen, 1978) and also observed by ourselves with a 


preliminary aluminum probe. 


32.5 bumped circuit” duplexer 

In our apparatus, a single coil is used instead of 
A Crossed COT! “coniteuration. “Thus, “1t* is necessary to 
isolate the detector from the output of the rf gated 


amplifier. One method is to connect the NMR coil to the 
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rf phase-sensitive detector by means of a quarter wave- 
length transmission line, with crossed diodes connected 
across the» input™“to the detector,“ During™a pulse) the 
diodes are turned on and therefore present a low impedance. 
The quarter-wavelength transmission line is thus shorted 
so that its input impedance becomes infinite. In the 
absence of a pulse, the diodes present a high impedance so 
that the transmission line is open-circuited. It therefore 
presents a low impedance and thus allows the NMR signal 
to pass to the rf phase-sensitive detector. 

The circuit described by McLachlan (1980) shown in 
Fig. 3.2 replaces the quarter-wavelength transmission line 
by lumped circuit components, and has superior performance 
to a transmission line. When a rf pulse is applied, the 
crossed diode pairs D,, Do and Dy conduct and Lo, Cy and 
Ro form a parallel tuned resonant circuit which presents a 
high impedance at the resonance frequency. The voltage 
appearing across Dg which are fast silicon diodes is about 


Or6o"V RMS. -"The*capacitorr’e iy con yunet ron “withthe "502 


3 
input impedance of the rf phase-sensitive detector, 
attenuates the voltage further so that only about 0.06 V 
RMS reaches the detector. When there is no rf pulse, all 
the diodes are open-circuited, and Lo and Cy form a series 
tuned resonant circuit. This presents a low impedance at 
resonance and allows the NMR signal to pass to the rf 


phase-sensitive detector. 


The purpose of the crossed diodes Dy is to connect 
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the rf gated amplifier to the NMR coil during the pulses, 
and to isolate it when the rf pulses are absent. The 
series resonant circuit is tuned by feeding a signal at 
the operating frequency into the probe terminal. The 
amplitude should be small enough so as not to turn on the 
crossed diodes. The capacitor Cy is tuned for maximum rf 
Signal entering the rf phase-sensitive detector. To tune 
the parallel sresonant ‘cireuit, switches Sy and So are 


closed, and Cy is varied to obtain maximum input impedance 


at the probe terminal. 


3.1.4 R.F. phase-sensitive detector 

NMR signals are detected by a Matec 615 tuned rf 
phase-sensitive detector. A cw level control adjusts the 
amplitude of the reference signal entering the rf phase- 
sensitive detector. For the work described in this thesis, 
the control was set at maximum at which setting the cw 
reference slevelwindieator read 17%. The unit is provided 
with high impedance and low impedance inputs. The low 
(502) impedance input was used in order to match the 
impedance of the tapped NMR coil. The unit is also equipped 
with high-Q and low-Q rf amplifier stages. The high-Q 
amplifier stage was used in our work because of the better 
Signal to noise ratio it provides. (The low-Q amplifier 
has a shorter recovery time than the high-Q amplifier. 
However, we were more interested in a good signal to noise 


ratio than in a short recovery time.) The unit is equipped 
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with a-control which enables the bandwidth to be varied 
between 0:9 kHz and 2500 kHz. For the work on a liquid 
sample, the bandwidth setting was B70, corresponding to a 
bandwidth of roughly 20 kHz. For the work on the crystals, 
the bandwidth setting was C70, corresponding to a band- 


width of roughly 200 kHz. 


a b-5> Ruileecacedfanplifier 

The Matec 515 rf gated amplifier, when plugged into 
the Matec 5100 gating modulator, provides gated rf pulses 
with amaximum power of approximately 1 kW. The output 
impedance is 502. The pulses from the pulse generator 
are introduced through the "External Modulation" input. 
The tampliiierars tuned *byeadjzustimerthéee'Rerer Tuning" 
control for*maximum output’ when®a 50 ohm’resistor is 


connected to the amplifier output. 


3.1.6 Signal averager 

The signal averager system consists of a Nicolet 
1170 high-speed signal averager with a Nicolet 171/2 
Signal digitizer plug-in unit. The data acquisition time 
can be varied from 1 us to 99.9 s per point. A post- 
trigger*sdelay ranging from 1 ws to 99.9 s' can be: introduced. 
Total memory consists of 2048 channels. The capacity of 
each channel is 20 bits. 

The signal digitizer plug-in receives analog signals, 


digitizes them and transfers the information to the high- 
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speed signal averager. The analog-to-converter 
TeSOLUtLON LS 6,79. OmelZzeDits ] The L2 -bite resolution 
was used in the present work. For preliminary work on 
a liquid sample, a Fabritek 1062 signal averager was 
used. However, the slow acquisition rate of 200 us per 


point made it unsuitable for work on the solid samples. 


3.1.7 Pulse generator 

The pulse generator used in our system was designed 
and constructed by the electronic shop in the Department 
Ole Nysices. eA schematic-dlacram, 1s. shown in Fig. 373. 
A synchronous down-counter is driven by a clock signal 
which is obtained from the frequency synthesizer. The 
initial settings of the down-counter are determined by one 
of four sets of decade thumbwheel switches, Pl, P2, P3 
and P4. The first set consists of three decade switches 
which can be set between O and 999. These determine the 
length of the first pulse between 1 and 1000 clock 
periods. (The extra clock period is caused by a time 
gelay an, the circuitry.) eine second set consists of a 


set of decade switches which can be set between 102 and 


9.999999x10°. These determine the separation time between 


the first and second pulse. The third set determines the 


length of the second pulse between 1 and Toe clock periods. 


The fourth set determines the pulse sequence interval, or 
the time between the second pulse and the first pulse if 


the pulse sequence is repeated. It can be set between Tor 
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PiAsure sais 
Schematic diagram of the pulse generator. The arrows 


indicate the direction of information flow. 
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and 9.9x10° clock periods. Since the operating frequency 
is 14.5 MHz, the pulse lengths can be varied up to 69 us, 
and the pulse separation and pulse sequence interval can 
be varied up to 69 seconds. 

At the beginning of a cycle, data is loaded from Pl 
into the synchronous down-counter via the thumbwheel input 
data selector. When the counter reaches zero, data gate 
Dieter closedv and data gated PZ) 1s opened so that’ the data 
from the P2 thumbwheel switches are loaded into the 
counter. When the counter reaches zero, data gate P2 is 
closed and P3 is opened and the procedure continues until 
the counter has counted down to zero for P4. The logic 
Signals are fed to the output buffer. 

Individual outputs corresponding to Pl, P2, P3 and 
P4 are provided. For experiments where single pulses are 
required, Outpuu © or PS 1s used. “in another output, Pl 
and P3 are added and this provides the gating signal for 
a two-pulse sequence such as the 180°-1t-90° sequence used 
in our experiments. The output buffer also provides narrow 
DiULses, which, occur at tho vend of each of Bl P2, P3 and. P4:. 
These are used to reset a pair of flip-flops in the rf 
gated amplifier for the reason to be mentioned in section 
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3.2 Details of Sample 


12,0 Description Ol sthessodium nitrate crystal 


The sodium nitrate single crystal used in the work 

was the same one as used by Spencer and Hughes (1978). 

It was obtained from the Harshaw Chemical Company, Cleve- 
land, U.S.A. and a spectroscopic analysis provided by 
Harshaw indicated the following impurities: Al, 3 ppm; 
Carob ile «3. DDIM: "Cue Iepom. co4).< dappm. ithe crysta 
sn cylindrical form,:J59mm an diameter and 11 mm in 
length. The three-fold axis of symmetry is accurately 


parallel to the end faces of the cylinder. 


Biase Crystallography. of NaNO, 
The Nitta Ce Pleo. NaNO, PS. SHOW el jer o.4e,) Laas 
material crystallizes in the rhombohedral system (space 


63q) with two molecules per unit cell. There 


group R3c, D 
are nitrogen atoms at each corner of the cell and at 1/2, 
Lio j l/ 2; the sodiumeatome  are;situated ate 1/4.) 1/4). 1/4 
ang 73/4). 3/4, 33/4 as ¢snownp an Fic. 3.4. ° The nitrate 
groups are planar, each nitrogen atom being situated at 
the centroid of an equilateral triangle of oxygen atoms 
whose plane is perpendicular to the three-fold axis. The 
sides of the oxygen triangles lie parallel to the projec- 
tion of the cell edges on a plane perpendicular to the 


three-fold axis, the orientation of the triangles 


alternating by 180° along the three-fold axis. There are 
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Figure 3.4 


Unit cell of. sedium nitrate 
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two sodium atoms, situated at sites A and B say, per unit 
cel bcot NaNO. crystal. These sodium atoms have a similar 
environment possessing 3 symmetry. However, the environ- 
ments are in each case differently oriented in space. 
These sodium nuclei have different intrinsic quadrupole 
relaxation probabilities, Win Wi p> Won? and Wop (Hughes 
and Reed, 1971). However, the time-averaged electric 
field gradient V (= Hoy = ~2V ), which determines the 
ZZ, 3% yy 

quadrupole splitting (Cohen and Reif, 1957), has the 
same value at the A and B sites. The resonances from 


23Na nuclei at these sites therefore overlap. Since 


mM=e3 pastor Sal the NMR spectrum consists of three lines, 
the so-called centre line and satellites (see section 2.2). 
Even though the A and B nuclei in general have different 
intrinsic quadrupole relaxation probabilities, they are 
strongly coupled by spin exchanges (Abragam and Proctor, 
1958; Bloembergen et al., 1959; Andrew and Swanson, 1960), 
which result from the overlapping of the resonances. Thus, 
the 23Na nuclei in NaNO. behave as a system of identical 
nuclei, if crystal imperfections are ignored. The 
effective quadrupole relaxation probabilities Wy and W 


2 


are the arithmetic means of Win and W and Woa and W 


Loa 2B 


respectively (Hughes and Reed, 1971). 


Syn (Preliminary? checksor athe tequipnent 


In nuclear magnetic resonance, the higher the 


frequency, the better the signal to noise ratio. Also, in 
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our work it was necessary to achieve good overlapping of 
the component lines of the quadrupole split spectrum. To 
achieve this, the second order quadrupole splitting must be 
small and this, in turn, means that the NMR frequency must 
be high (see Eqs. (2.43) and (2.44)). However, the maximum 
usable frequency is determined by the magnetic field 
strength available. We worked with a field of approxi- 
mately 13 kilogauss, which is close to the maximum field 

of our magnet. The gyromagnetic ratio of the ae NE nucleus 
is such that the NMR frequency in a field of 13 kilogauss 
is 14.5 MHz. This is the frequency for which the sample 
probe was designed (see section 3.1.2). 

The performance of the pulsed NMR spectrometer was 
checked using the CoA Signal from a 2 ml sample of a 
doped saturated aqueous solution of sodium nitrate. The 
doping consisted of approximately M/10 concentration of 
paramagnetic Fe(NO,). in order to reduce the spin-lattice 
relaxation time to the order of 10 ms. This is done so 
that the spin system could recover quickly from the effects 
of a pulse. This allows pulsed experiments to be repeated 
at frequent intervals. 

The gain of the rf gated amplifier was turned up to 
maximum, at which setting the voltage at the input to the 
NMR coil (across the single turn) was 600 volts peak-to- 
peak Os This corresponds to approximately 10 KV Sp 


across the whole NMR coil. All the measurements reported 


in this thesis were carried out using this maximum pulse 
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amplitude. The FID immediately following each pulse was 
observed on an oscilloscope connected to the video output 
of the rf phase-sensitive detector. The oscilloscope was 
triggered by the trailing edge of the pulse. The repeti- 
tion rate of the pulses was 1.4 seconds leaving plenty of 
time for the spin system to relax back to equilibrium 
between pulses. The strength of the magnetic field was 
varied until a proper FID was observed; if the field is 
off-resonance, the nuclear precession frequency is 
different from the frequency of the reference signal fed 
to the rf phase-sensitive detector and the observed signal 
has the torm of a damped oscillation. 

After the proper FID was obtained, the reference 
phase was adjusted in order to obtain a maximum signal. 
The length of the rf pulses was then adjusted until the 
amplitude of the FID was a maximum, i.e. until the pulse 
was a 90° pulse. The pulse length indicated by the 
setting on the pulse generator was approximately 115 rf 
periods or 7.97 Us.” "Phis pulse length is*quite short, 
bearing in mind the rf power available. It shows the 
effectiveness of our tapped-coil impedance match. 

Careful observation of the amplitude of the FID 
following a pulse of approximately 90° showed an inconsis- 
tency. At certain pulse lengths, every fourth pulse was 
found to be four rf periods shorter than the others. The 
problem was overcome by resetting a pair of flip-flops in 


the rf gated amplifier to coincide with the beginning of 
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the pulses from the pulse generator, as suggested by Mr. 
KE. Cairns and Mr. W. Siewert. This required a small 
modification to the circuit of the rf gated amplifier, 
and the provision of narrow 'reset pulses' obtained from 
the pulse generator. 

In order to obtain a 180°-tT-90° pulse sequence for 
measuring spin-lattice relaxation times, the length of 
the first pulse generated by the pulse generator was 
increased from zero until the amplitude of the FID 
immediately following the pulse was zero. The absence of 
a FID means that there is no transverse magnetization and 
tris implies that the pulse 1s 2° 180 "pulse. The Teneth 
of the second pulse generated by the pulse generator was 
increased from zero until the amplitude of the FID 
immediately following the pulse was a maximum, i.e. so 
that the second’ pulse is a 9O°~ pulse.” The’ length of the 
180° pulse was found to be 205 rf periods or 14.2 us. 
This is’ sliehtly less than twice» the length of the 90° 
pulse. This is presumably because of phase shifts which 
occur at the leading and trailing edges of each pulse. 

The operation of the equipment was checked by 


a EA nuclei in the doped saturated 


measuring Ty of the 
aqueous solution of NaNO.. These measurements were made 
using the Fabritek 1062 Signal Averager, operating at its 
Tastes: data acgauisition rate of “200 s/point. 


As described in seccion 2.1.2, the approach to 


equilibrium of the z-component of magnetization following 
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a 160 "pulse will be of the’ form 


Qn (Mo-M,) =A - (t/T,). C3 2tl)) 


Experimentally My can be found by measuring the amplitude 
of the FID following a 90° pulse and M, can be found by 
measuring the amplitude of the FID following a 90° pulse 
Anpried at a vime tT @iter 2 160. pulse. “Since” the 
difference between My and M, is required, it 1s not 
necessary to measure the baseline of the FID, provided the 
baseline is the same in both cases. Examination of the 
baseline showed it to be quite stable, so that it was 

only necessary to measure the difference between the 
magnitude of the FID's associated with My and M.. These 
values were measured in the fifth channel of the Fabritek 
Signal Averager, i.e. 5 times 200 us or 1 ms after the 

Jo "pulse. “LoO™optaim cOod sienal to Noise ratio, signals 
were recorded 32 times in the signal averager. In order 

to minimize any error due to a change in gain of the system, 
values of My were measured between each M, measurement. 
Thus, the sequence was Mo» M(t); Mo Mo(T)); 

M(t )> Mo: The arithmetic mean of the two Mo values on 
either side of a MS value was used in determining &n (Mo-M,) 
fr ea. Coby therdara enown i Cis, o.0 talbis On 2 
straight line, in agreement with Eq. (3.1). The value of 
Te-Tound £rom 2° lLedst suuares tit Of the data toa Straight 
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line is 20.4 ms. Repeating the measurements gave values 
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ranging from 20.0 ms to 20.8 ms, showing reasonably 
Satisfactory reproducibility. 

According to Andrew (1955), the spin-lattice 
relaxation time of protons in a M/10 aqueous solution of 
Perr ice Ons SY Lo, ise) eoance T) is inversely proportional 
to the square of the gyromagnetic ratio of the nucleus 
involved, and since the gyromagnetic ratio of the proton 


1s-3./8 times as large as that, of the ona nucleus, the 


expected Ty of 439 in a M/10 aqueous solution of ferric 
go0nsS ws 16.5 ms. Bearing in mind’ the uncertainty in the 
concentration and the simplifications made in the theory, 
the agreement between experiment and theory is satis- 
LAcCtOry. 

Examination of the My values obtained with the 
doped aqueous solution of NaNO, Showed no evidence of any 
systematic variation of My With “time. “Thus, 10 was 
decided in future work to use the overall mean value of 
Mo» rather than the mean of the Mo values measured 
immediately before and after each M.. 

The data shown in Fig. 3.5 were taken with tT 
increasing monotonically throughout the experiment. This 
Le not Food practice, and it was decided in future to 
choose the tT values in a random order during the experi- 
ment. In this way, systematic errors due to a time-varying 
gain or baseline will show up as an enhanced scatter of 


the data points, rather than as a curvature of the data 


points in the semilogarithmic plot. 
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Figure’ S75 
3 , a3 f 
ee versus pulse separation time Tt for Na in 2 


doped saturated aqueous solution of NaNO. at 299K. 
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Examination of the data also showed the error in 
M, was due to noise and was independent of the value of 
t (and was also the same.as the error in M,): Thus, the 
error in the ee is the same for each data point. 
However, the error in &n(M.-M,) is inversely proportional 


to (M-M,) as can be seen by differentiation: 
6 &n(M-M,) = 5(M\-M,)/(M-M,) - C32) 


Since (M-M,) varies by about an order of magnitude ina 


typical T, experiment, the random error associated with 


i 
each data point also varies by about an order of magni- 
tude. It is therefore important to perform a weighted 
teast squares fit,to, the datas It is also helpful to 
determine the error in T) obtained from the least squares 


fit. The remaining T+ values to. be reported in, this thesis 


1 
were therefore obtained by means of a weighted least 
squares fit of the data. Also, two error estimates were 
obtained. One is derived from the scatter of the data 
about the straight line. The other is derived from the 
standard error associated with each data point. An 
estimate of this standard error can be obtained from the 
scatter in the various Ma values. The formulae required to 
perform a weighted least squares fit and to calculate 

the errors in the fitted parameters, are presented in 


Appendix 1. The method of obtaining the standard devia- 


tions and weighting factors needed in our least squares 
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fitting procedure are described in Appendix 2. A program 
for the Hewlett-Packard 67 programmable calculator that 


would calculate T; and. the error in TS from values of 


(M-M,) and t is described in Appendix 3. 


3.4 Measurement of the spin-lattice relaxation time of 


aera in sodium chloride 


In order to compare our T) values with those of 


other workers, and to look for possible systematic errors, 


we measured the spin-lattice relaxation time of 23Na at 


a Single crystal of NaCl at room temperature. Even 


though this T, is believed to be independent of crystal 


iL 
orientation relative to the external magnetic field (Zak, 
1964; Snyder and Hughes, 1971), we aligned the crystal so 
that H was along the 100 direction. 

Since the T, value is approximately 14 s (Schumacher, 
1958; Satoh, 1965; Spencer and Hughes, 1978), the pulse 


sequence interval was set at its maximum value of 69 seconds 


Or avouwr OL so that the starting magnetization is 


17 
(1-227?) Or 0, 96s Mo: Because of the long time between 
pulse sequences, the signals were averaged only 8 times, 
even though this sacrificed the signal to noise ratio. 
According to the Matec manual, our rf phase-sensitive 
detector should be linear provided the video putput is less 
than 1 Lae We kept well below this limit by adjusting 


the receiver gain so that the maximum output was only 0.25 


re The bandwidth in the rf phase-sensitive detector was 
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Set. at the position C70 corresponding to a bandwidth of 
Frougchty 200 kHz. 

Since the decay of FID of 23Na in NaCl at room 
temperature is much faster than the decay of the FID of 


23 


Na nuclei in the doped aqueous solution of NaNO a 


3 
much faster data acquisition rate than the 200 us per 
point of the Fabritek signal averager was needed. We 
therefore used the Nicolet 1170 signal averager, borrowed 
Brom the University of Caleary. dor this work. A data 
acquisition rate of 1 us per point was used. Also, a 
post trigger time delay of 28 us was introduced so that 
the rf phase-sensitive detector output during the 
recovery period was not recorded. The t values were 
selected in random order for the reason discussed in 
section 3.3. The FID signals were recorded in a quarter 
of the memory of the signal averager or 512 channels. 

The data shown in Fig. 3.6 was obtained from the 
FID's recorded in the seventeenth channel of the signal 
averager. Thus, they were recorded 17+28 = 35 us after 
the trailing edge of the 90~- pulse. For an ideal crystal, 
the shape of the FID should be independent of t. Thus, 
it should not matter where (i.e. at what time) the magni- 
tude of the FID is measured, provided it is done 
consistently for all values of 2. (Ciraditionally; the 
magnitude of the FID has been measured as soon as possible 
after the 90° pulse). However, in order to check this 


point, we measured the FID's in the fourth, fifty-seventh 
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Figure 3.6 
Mo-M, versus pulse separation time t for oe in NaCl 
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and ninety-seventh channels, in addition to the 
seventeenth channel. The values of T; were also obtained 
from the same FID's using data in the fourth, fifty- 
seventh and ninety-seventh channels. (These channels were 
far enough apart so that little correlation between the 
values is expected). The calculated T, values are shown 
in Table 1, together with the standard (68% confidence) 
error associated with the scatter and the standard error 
associated with the standard deviation per point. These 
two types of error estimates have been called external 
and internal errors respectively (Birge, 1932). It can be 
seen from the table that the two types of errors are of 
Similar magnitude, indicating that the proper fitting 
ETNictiOniisn a) Strapehhdanes eliorapurpose of calculation, 
the larger of the two errors was used in each case, as 
shown in the fifth column of Table 1. The weighted mean 
TS is 14.6+0.15 seconds. This is compared with other 
published values in section 4. 

Because the magnet enclosure was slightly warmer 
than room temperature, a thermometer was inserted into the 
sample probe (with rf pulses turned off) and the tempera- 


ture was found to be 26°C: “This was taken to be ‘the 


sample temperature. 
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Table l 


, ona Sin eNaGl woot ied usineedata ftom 


different channels of the signal averager. 


Channel number Ty Error Error T 


(Seatter)ice(Stdandev.)soGseconds) 


4 14.37 0.20 Or22 TA S720 22 

tar 14.91 OF 30 0-20 149 FPO 230 

auf 14.76 0.34 0.30 TAO ros 

96 eos Oak 0.42 LAS OteOl 42 
Weighted mean T, = 14.60+0.15 seconds. 
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3.5 Measurement of the orientation dependence of the 


spin-lattice relaxation time of oon in NaNO... 


S20. Installation. of the NaNO, crystal in the sample 
holder. 

A special sample holder was constructed in such a 
Wayerllat tCOULd De 'Si1d into (Ghe NMR. coal and: ‘so that 
the crystal could be rotated independently about the two 
axecs-shown-in Figs 3.7. Axisel is the axis of the sample 
holder and the NMR coil. It: as perpendicular to He 
Axis 2 is perpendicular to axis 1. The crystal A was 
mounted in a collar B which can be rotated about axis 2. 
The collar is equipped with a vernier scale (see Fig. 3.7) 
so that its orientation can be read with an accuracy of 
about one degree. The sample holder can be rotated about 
axis 1 in the sample probe. The probe is equipped with 
a circular scale C and a vernier which can be read to 
within one tenth of a degree. 

The crystal was aligned in the collar so that the 
three-fold axis was along axis 2. This was done in the 
following way. It can be shown from Eq. (2.44) that the 
separation of the two oA satellite lines in NaNO, ee 
a maximum at 6 = O (and 180°). The crystal was set in the 
collar so that the end faces are perpendicular to axis l. 
This 18 the situation suown in fie. 3.¢. The: sample 
holder was then rotated about axis 1 until the satellite 


separation was a maximum, i.e. until 6 was 0°. This was 
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achieved by setting the magnetic field to be on resonance 
for one of the satellites and rotating the crystal about 
axis 1 until the resonant magnetic field was a maximum 

or minimum, depending upon the satellite chosen. The 
erysital Gand the collar), was then rotated, through .90° 


about axis 2 and the sample holder was again rotated 


about axis 1 until the satellite separation was a maximum. 


This satellite separation was less than the value corres- 
ponding to 8 = O° because axis 2 and the symmetry axis 
were not parallel to one another. The crystal was then 
rotated in the collar and the above measurements were 
carried out again. The process was repeated until there 
was no difference in the satellite separation when the 
Grvstal wasuroteted throughs [On sanout. axis 2 sAX1 Ss) 2 and 
the symmetry axis were then parallel to each other. 

An important consequence of this adjustment is that 
rotation about axis 1 changes 6 but not 9¢, whereas rota- 
Mion 1bOUt ax LS) .2 sCOanC eS. 40 DUT ANOG Os. sshnisSesi tation 
greatly facilitates the measurements to be described in 


the remainder ot this ~thesis. 


3.5.2 Choice of crystal coordinate system for NaNO... 

In choosing the crystal coordinate system referred 
tO in esection 2.3 s0-asito makewtiull use of crystal 
symmetry, the z axis is selected to coincide with the 
three-fold symmetry axis of the crystal. Also, the x 


axis is selected so that it lies in a glide plane. 
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(Each glide plane contains the three-fold axis and a 
uid. te cert” Cage). #" Thee location of tne xX axis anethis 
crystal had previously been determined by means of an X-ray 
precession photograph of the crystal (Spencer and Hughes, 
1978). However, there is an ambiguity of 60° in the 
orientation of the x axis which is not resolved by the 
X-ray method. This ambiguity is removed by arbitrary 
selecting the (right handed) coordinate system xyz so that 
the unit cell edge which lies in the xz-plane extends from 
the origin into the second or fourth quadrants (rather than 
inca, thesdirst and third quadrants). The orientation. of 
this axis was found using NMR, by comparing the width of 
the SON resonance’ lime” at the-setting @°="54.7°° and 
o = O° and 60°. The second moment of the centre line in 
the frequency domain, due to the magnetic dipole-dipole 
interaction, -1s*given by “CAndrew*vet al. 7" 1962)" 
S = {29.57-2.97c0s20-9.14c0s40+6.04(2sin26-sin46)cos3$}10* Hz" 
Gaia) 

relative to the crystal coordinate system defined above. 
At 6 = 54.7°, the satellites overlap the centre line. 
However, the ¢ dependence of the second moment of the 
combined resonance line should still be of the same form as 
Bay, (360). 

It can be seen from this equation that the dipolar 
second moment for 0 = 54:7° is a maximum at 9 = 0° and a 


minimum at ¢ = 60°. The FID which is the Fourier Transform 
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of the NMR signal in the frequency domain, should there- 
fore be narrower in the time domain at 6 = 0° (and 6 = 
54.7°) than at @ = 60°. The “°Na FID's in NaNO, obtained 
fOr. 9 —2354 7 and o = 0. sand ve0- are ‘shownean Hie. 3-3. 
The FID's were averaged eight times and then read out 
from the signal averager via the X-Y recorder. It is clear 
from the above discussion that the @ values associated 
Witne the two FID"s are as shown in Fig. 3.6. This 
uniquely identifies the orientation of the x axis in our 
Crys ual. 

3.5.3 Measurement of T) of oN in NaNO... 

The NaNO. crystal was oriented so that o was equal 
EOvoO VOY Totatins bhe coOllareby ADOUt! AxTS (2GSseer ilo. 6357.) + 
The angle 98 was set at 54.7°, where resonance overlap 
occurs, by rotating the sample holder about axis 1. The 
position of best overlap was taken to be the-orientation 
at which the FID was approximately Gaussian in shape. 
Rotation of the crystal through) about 021° inveither 
direction caused the FID to change shape significantly and 
in particular to overshoot because the satellites became 
off resonance. Thus, the setting of the crystal was very 
delicate. 

The equipment was allowed to warm up for one hour 
before Ty data were taken. Since the T) of oN nuclei is 


approximately 6 seconds (Kasahara et al., 1977; Spencer 


and Hughes, 1978), a repetition rate of 34.5 seconds was 
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Figure 3.8 


Free induction decay F(t) of 23a in NaNO. ate @ = 54572 
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used. Also, eleven different t values ranging from 0.1 to 
15 seconds were used in random order. Otherwise the method 
was the same as for the NaCl crystal. In particular, each 
FID was averaged eight times. Data obtained in the fourth 
channel (32 us after the end of the 90° pulse) are shown 
Tie h1e7go.9. 8 As was mhetcase tior NaCl. T, values were 
also obtained from the data in the seventeenth, fifty- 
seventh and ninety-seventh channels, and the weighted mean 
was calculated. The measurements were repeated, using the 
same tT values, with » covering the range 0° to 120° in 15° 
steps. For reasons discussed in section 4, great care was 
taken to keep the operating conditions, and in particular 
the control settings on the rf phase-sensitive detector, 
the same for all 9 values. The Ty) values are summarized 
iieebabile 2: 


The temperature of the sample was 299K, as in the 


case of the NaCl crystal. 
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Figure 3.9 


M-M, versus pulse separation time t for 23ye in NaNO, 


at 299K and sor =) 3o0re 
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CHAPTER 4 
DISCUSSION 


We found that Ty) of 239 in Our Single Clyseal OL 
NaCl at room temperature (299K) was 14.60+0.15 seconds, 
where the error limits include random errors but not 
necessarily all types of systematic errors. Other 
published values of Ty) for oon, an NaCl are 14).220.5 
seconds at 300K (Schumacher, 1958), 14.2+0.5 seconds at 
SOU Coacoh, 19600), and 14. 12074 Seconds at 290K 
(Spencer and Hughes, 1978). (Satoh's value has been read 
Gita graph. He’ did not explLicitiy state his error: Jimits:. 
However, it can be inferred that they were approximately 
+0.5 seconds.) The TS of NaCl is approximately inversely 
proportional to the square of the absolute temperature. 
Spencer and Hughes' value therefore corresponds to 13.75+ 
0.2 seconds at 299K. While our value agrees with values 
found by Schumacher and Satoh within the combined error 
limits, it does not agree with the value found by Spencer 
and Hughes (1978) using a cw technique. Satoh, and 
Spencer and Hughes, used a single crystal of NaCl, as we 
did in our measurements, so the discrepancy cannot be 
attributed to using a different type of sample. The reason 
for the discrepancy between our value and that of Spencer 
and Hughes is not known. However, the fact that our T, iS 
longer than all the other published values raises the 
question whether there is a systematic error in our TS 


values, due to nonlinearity of the rf phase-sensitive 
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detector for example. 


The T) of a in NaNO. 


different channels (at different times after the 90° pulse), 


at 299K, as measured in 


arelshown?t formvargous.outoumne Tablengre Vatues oF Ty» 
calculated from data obtained using different channels, 
are in quite good agreement with each other, except for 

@ = 60° where there is a significant tendency for T, to 
increase with increasing channel number. 

The resonance linewidth is a minimum at @o = 60° 
VAndrewvet ali, 1962)" (see also section 3.5.2). The 
resonance overlap achieved will therefore be less at this 
orientation than at any other. However, if this is the 


reason for the dependence of T, on channel number, one 


1 
would expect almost the same trend at 0 = 45° and 75°, 
and there is little evidence of this in Table 2. 


Another possible reason why the measured T, might 


1 
depend upon the channel number is the presence of crystal 

imperfections. In an imperfect region of the crystal, the 
satellites will not overlap the centre line properly, and 

will not “be scovtpléedsevobit by spin*exchanges? “Inthatycase , 
the centre line and satellites will relax independently of 
one another. This is shown mathematically by the following. 


In the case of a common spin temperature, the relaxation 


of the population differences N, between the m = 3/2 and 


1 
1/2 levels, No between the m = 1/2 and -1/2 levels and 
ry between the m = -1/2 and -3/2 levels following a 180° 
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where ny is the population difference between adjacent 
levels in thermal equilibrium. If there are no spin 
exchanges between the centre line and satellites, the 
relaxation following a 180° pulse is given by (cf. Andrew 


and mMunstal 17u961) 
N= n [1-2 expl(-2 Wot) ] (4,2) 


N => he (3 +2rme xpi Wet a P40 expaié -2W ate 
Oo O 1 2 
(4.3) 
In imperfect regions of the crystal, the satellites will 
decay more rapidly than the centre line in the FID 


foliowdne ai 90°d pulser? Thusye thei isatellites!) will? contri= 


bute less to the FID as time increases, i.e., as the 
channel number increases. Values of Ty) would then in 
general depend upon the channel number, though if W, and 


a 


Wo are’ equal ,Jat canvbe*seent? from Eqs. (4.1)%°(4. 2)" and 


(4.3) that the relaxation is the same whether or not spin 
exchanges occur. It can be shown, using the M-tensor 
components found by Spencer and Hughes (1978), that at 


9 23 


23) 
$ase 0leanndii1217 8! abe de S2o0ee Imperfections should thus 


have the largest effect on the measured T, values at ¢$ = 


1 


60°. However, it has been shown (Hughes, 1981) that even 


if imperfections were such that the centre line alone was 


bAeice Wo/W, for Na in NaNO, varies between 0.967 at 
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observable in 10% of the volume of the crystal, the 
effect would be to reduce the T; measured at $9 = 60° by 
less than 1%. A rough comparison of the amplitude of the 
FID ot the centre bine only (Cat © = 0°) with that of the 
combined satellites and centre line (at 8 = 54.7°) showed 
that the volume of the crystal affected by imperfection 


is small. Thus, the dependence of T, upon channel number 


1 
at ¢ = 60° is unlikely to be due to crystal imperfections, 
and the reason for it is unknown. 


The weighted mean values of T, given in the last 


if 
column Of Table 2° are “shown/ac-a function of “6. in Figw 4.12 
The data look as if they vary as -sin3o. However, 

thes is probably coincidental: > When the error bars’ are 
taken into account, the data fit quite well to the dashed 
straight line drawn at the weighted mean value of T = 6.19 
seconds. Of the nine’ data points; only four ditier ‘from 
6.19 seconds by much more than the error bar. Bearing in 
mind that the error bars represent 68% confidence limits, 
that is not too unreasonable, especially since it has been 
found in many laboratories that the day-to-day reproduci- 


bility of T, values measured by the pulse technique is not 


1 
high. We therefore deduce that T) for AS in NaNO, is 


independent of » as predicted theoretically (Hughes and 
Spencer, 1979) (see also section 2.3). As will be 
mentioned later, there is evidence of a systematic error in 
our values of T,. However, the experimental conditions 


1 


were kept the same for all values of ¢. The systematic 
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Figure 4.1 
T; of oye in. a single crystal o: NaNO. at 299K, as a 
function of > for 8 = 54.7°. The error bars represent 


the larger of the two standard errors calculated by the 


least squares fit. 
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error should therefore Mot rafifect the iformyvof the 


orientation dependence of Ti: 
Our TS value of 6.19 seconds for 23y8 in NaNO. 


299K can be compared with two other values. Using fitted 


at 


parameters found by Kasahara et al. (1977), we find a T, 


value tat e299K eofe6 OSuasecondsslvBy substituting yalues 


of the M-tensor components found by Spencer and Hughes 


(1978) in Eq. (2.57), we find that the T, of BOR ais 


NaNOF ato vbed6Gt42+000S8 secondsnat (295K sa0Correcting this 


a 


value to 299K using the temperature dependence of T; 


Poundeby *Kasaharacerial va i977) towe: find that T, at 299K 
should be 6.2340.08 seconds. Our value of 6.19 seconds is 
in excellent agreement with the values deduced from Spencer 
and Hughes' work and from the measurements of Kasahara 
et al. There is no evidence of the kind of systematic 
discrepancy with T) of eon pntNaCitntHowevers astwitl 
become clear in the next paragraph, this excellent agree- 
ment is probably coincidental. 

Finally, we consider possible defects in our equip- 
ment or procedure that might cause erroneous values of T,- 
The measurements reported in this thesis were obtained with 
a smaller reference signal (17% on the cw reference level 
meter) than that recommended in the manual of the rf phase- 
sensitive detector (50% on the cw reference level meter). 
This was because of a Joss of signal in the phase shitter. 


Recently, measurements made by Pandey (1981) indicate that 


the use of a reference signal of the recommended amplitude 
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gives T) values which are roughly 2% larger than those 
reported in this thesis. Therefore, this cannot explain 
the fact that our measured value of T, for Te) in NaCl 


i 
was longer than other published values, Pandey (1981) 


also found that use of the low-Q setting in the rf phase- 


sensitive detector, gave T, values approximately 3% 


ug 
larger than the values obtained with the high-Q setting 

used for the measurements reported in this thesis. 
Generally, quite good agreement is found between T, values 
measured on different days with the same control settings 

on the rf phase-sensitive detector. All this suggests 

that the rf phase-sensitive detector may not be perfectly 
junear-oand, thatdtheinonlineari ty? varies with: the ‘control 
settings of the rf phase-sensitive detector. Some confirma- 
234 


tion of this has been found by reevaluating the T, On 


usineytthe .asineix data. points! (those with +t 2.5 


a 


in NaNO. 


seconds, see Fig. 3.9). These T, values were 1% smaller 


on the average than those found by taking all eleven data 
points. This suggests that the equipment may indeed be 
nonlinear. In order to achieve good accuracy of say +1% 

in TS values, our rf phase-sensitive detector may therefore 


need to be replaced by a better instrument. 
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APPENDIX 1 


Weighted least squares fit to a straight line, and the 


error in the fitted parameters. 


We suppose that n data points X1Vy> XoVor ces sh Pi 
are to be fitted to the straight line 
Vos iter cs CALs) 


We suppose that the individual points have associated 


(relative) weighting factors w ees According 


sip: Wino 
to Squires (1976), the best value of the slope m is 


given by 


ae) w,(*,-X) Oe @abeey 


where 


2 -,2 
ARES iat OO) (eAl,.3) 
Since the mean value x is only available after all the 
data have been inserted and the fit has been completed, 
it is more convenient to express x in terms of 


summations as follows: 


x =) wX,/) Wi: GAleea) 
it -Lollows. that 


Dive y wox,7 Set) w x.) 7/) W,- GATS) 


Substituting in Ed. (Al v2), we tina that 
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This expression agrees with Eq. (6-12) given by 
Bevington (1969), provided Wi is replaced by 1S 
where Oo, is the’ standard error of the ith data point. 


Similarly, we find that the best value of c is 


given by 
Cc = y - mx (A1.7) 


or 


= DMN GE 4 


2 
: Lw iV LW 4%5 


5 5 CALS) 
iwydrye, — (hw i%)) 
The standard error inm and ec, as determined 


from the standard error of each data point, is given 


by (Bevington, 1969) 
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ios blow (xe) = Clay ow 
(A1.9) 
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a. 
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Ge fg deus Guicaye ert sad cyh fds 
(Al. 10) 


Q 
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An equally good, or perhaps better, estimate of 
the error is provided by the scatter of the data about 
the straight line. The standard error in m and, as 


determined by the scatter, is given by (Squires, 1976) 
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Z -,2,% 
= d 2 ca 2 
Ol Cowes mena 2) Lw (x, x) ey CALAI1) 
and 
a. = {[Cl/lw,) + x-/)wicx.-) lw.d.7/(n-2)}2,  (AT.12) 
G ah yea. a: : . 
where 
d. senna Tk eace CAT. Is) 
mS the deviation Of the ath data point Cin the y 
Gurection)) from the straieht bine. Substituting 
for x and d. usin Bas. WCAl4) sand CAL 13): we can 
express a and ao in terms of summations as follows: 
as 
tw. (w.y, eo }w, -m“Jw.x, “-2e)w.y.) ; 
= ligase WE i ish rea 
Naa Ent ee SOR ns Oe OG ae SO ee el GNIS) 
(n-2)(}w, )w x, = ()w5x,) ) 
and 
2 pe Zi 2 z 
A: ()w5x, eT aye ee) ee) ee -2c)w,y,) 
c 2 2 : 
(n-2) [lwidwix,  - (lwyx07] (A1.15) 


We have not found these expressions (Al.14) and (A1.15) 
in the literature. However, they are very useful for 
calculating the standard error in the fitted parameter 
m and c by means of a programmable calculator, since 


the data need only be retained in the form of summations. 
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(The only disadvantage oieueinc bas. CAl.44) and, CAI. 15) 
is that large numbers of nearly equal magnitude are 
subtracted. Care must therefore be taken to use enough 
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APPENDIX 2 


Method of obtaining the standard deviations and weighting 
factors needed in our least squares fitting procedure. 


Experimentally we found that the standard deviation 
caused by random error in MS. is independent of Tt and is 
therefore the same as that associated with Mo: This 
standard deviation, o[M.], was calculated from the 
scatter of the values of Mo: Also, the mean value of 
Mo» say Mo was used in calculating values of &n(M-M,). 


The standard deviation of &n(M-M_) is given by 
M - = vem Na <2 
o[£n(M. M1 oLM. Mo1/(M, M,)- Az 1) 


If we assume that the errors in M, and M, are uncorrelated, 


then the standard deviation of &n(M-M,) will be given by 
ofan(M-M_)] = {Coli.])? + Colm_])7}2/(i-M_) (A2.2) 
en pie O zy Oy was : 


The standard deviation in M,, viz. o[M_], is equal to the 
o[M.] obtained previously. The standard’ deviation of My» 
viz. o[M] is given by oLM J/yn, where n is the number 


ofeimdividual values of Mo used to calculate Mo: (For 
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our measurements of T, of Na, in a single crystal of 
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NaNO., n wasre2s jewilhus.- whe.cstandard deviaiaoman 


&n(M.-M,) is 


o[ 2n(M.-M,) J = o[M.] he ely ney (M.-M,) - CAO. 3) 
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The weighting factor associated with each data point is 
inversely proportional to the square of the standard 
deviation. Thus, the weighting factor associated with 


the ith data point is given by 
w, © (M -M_)7/(olM_])? {1 + (1/n)} (A2.4) 
i Cuens fe) ‘ 


We obtain the value of T, by performing a least 


squares fit of &n(M_-M,) versus T, and finding T from the 
Slopes Odesther best «straieht dine wsing ig <:CAle6) and The 
ren i ee OY Oma a T) isy found! usdng Haste:(Ade9o) and (Aled4)£ We 
manually calculate the standard deviation of (M_-M,) which 
is given by the numerator on the right hand side of 
Bote GAZ.) .teelhas.ouantity, which is the sameefor all 
data points, is then fed to a programmable calculator 
which performs the least squares fit. The program 
evaluates o[&’n(M.-M,)] for each data point by means of 
Eq. (A2.3). It then takes the inverse square of these 
quantities and uses them as the weighting factors in 
evaluating the summations which occur in Eqs. (Al1.6), 
(Al.9) and (Al1.14). The standard deviation Oo. which 
appears.ine Ede Al«,.9)4 asdofycourse o[&£n(M-M,) I], and this 
has already been calculated. 

The program used to. calculate T. is described in 


sl 
Appendix 3. 
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APPENDIX 3 


Program Lor calculating iy -trom values’ o. (Mo-M,) Wale Pandy 


A 


It can be’Seen from Eqs.%(2+19) and (Al.1) that 


Tor=° -l/m. ALSOP On the standard error in T,> is given 
dk 


Dy O = 0 ine Two estimates of a are obtained, 
ie oo eal 


corresponding to the two different estimates of On given 
by Eas. (Al .9)@amd@ (Al .14). For conveniencergdn plotting 
mhe data craphically; the progamjpaiso computeszc, the 
intercept of the straight line, and the corresponding 
value of (Mo-M,). It also caleilates the value of +t for 
which (M,-M,) msureduced by exactly a factor*oi 10. 

The following summations are accumulated in the 


registers indicated: 


lw, Reg. 0 
lw 5X, Reg. 1 
Jw x7 Reg .2 2 
1s ie (where yes &n(M-M,) 4) Reg. 3 
ye Reg. 4 
lwiy.- Reg. 5 


yg 
The standard deviation o[M.-M, 1], which has previously 


been manually determined, is stored in Reg. 9. 


The program, which was developed by my supervisor 


fom the HP 67 programmable calculator, is asHollows: 
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To obtain the value of Ty > the programas hirst 
loaded. Then the A key is pressed. This clears the 
registers in preparation for the calculation and sets the 
Conurect display. «At this point, the value of o[M.-M_] is 
inserted in Reg. 9. The first value of Tt is entered, 
followed by the ENTER key. This is followed by the first 
value of Me The B key is then pressed and this 
evaluates the summations associated with the first data 
point. The second value of tT is then entered, followed 
by the ENTER key, the second value of (M\-M,)> and the B 
key in that order. The process is repeated for all the 
other data, atowhich point; the summations in Rees. 0-to 5 


are complete. Key C is pressed and this causes the value 


of T, to be computed and presented in Reg. A. The external 


1 
error 1 TS CeLTroOreaie to ther secattervot (the dara, points) 
isepresentedein Ree. B. The internal.error in T.; Cerror 


1 


associated with the standard deviation of each value of 
in(M -M_)) Jsepivensinenec ©. “The best value of MM, at 
t = O is given in Reg. D, and the value of T at which 


(Mo-M,) has fallen by a. factor of ten is given in Reg. E. 
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